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We present results of lattice dynamics calculations of Poisson’s ratio (PR) for solid hydrogen and 
rare gas solids (He, Ne, Ar, Kr and Xe) under pressure. Using two complementary approaches - the 
semi-empirical many-body calculations and the first-principle density-functional theory calculations 
we found three different types of pressure dependencies of PR. While for solid helium PR monoton- 
ically decreases with rising pressure, for Ar, Kr, and Xe it monotonically increases with pressure. 

For solid hydrogen and Ne the pressure dependencies of PR are non-monotonic displaying rather 
deep minimums. The role of the intermolecular potentials in this diversity of patterns is discussed. 
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At low temperatures and pressures solid helium is an 
ultimate quantum solid displaying such phenomena as 
zero-temperature quantum melting and quantum diffu¬ 
sion. As atomic masses and interatomic forces increase 
in the sequence Ne, Ar, Kr, and Xe quantum effects in 
their properties become progressively less pronounced. 
Solid hydrogen is the only molecular quantum crystal 
where both translational and rotational motions of the 
molecules are quantum. Translational quantum effects 
decrease with increasing pressure. 

Quantum and classical solids respond to the applied 
pressure differently. When pressure is applied to a clas¬ 
sical solid the atoms are ’’pushed into” the hard cores 
of the potential; as a result of this core, the compress¬ 
ibility is usually quite small. Typically, the pressure of 
1 GPa results in a few percent change in molar volume. 
At the same time, quantum solids hydrogen and helium 
are highly compressible. For hydrogen the pressure of 1 
GPa results in a 100% change in volume. The physical 
reason for this is that the lattice is highly blown up due 
to the zero-point kinetic energy. The initial compression 
works against the weaker ’’kinetic pressure” rather than 
the harder ” core pressure”. 

One of fundamental thermodynamic characteristics de¬ 
scribing behavior of a material under mechanical load is 
Poisson’s ratio [l|, Q. For isotropic elastic materials the 
Poisson’s ratio is uniquely determined by the ratio of the 
bulk modulus B to the shear modulus G, which relate to 
the change in size and shape respectively^: 
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As can be seen from this equation, PR can take values 
between -1 {B/G 0) and 1/2 {BjG cxo). The lower 
limit corresponds to the case where the material does not 
change its shape and upper limit corresponds to the case 


when the volume remains unchanged. Materials with 
small PR (small BjG), such as cork, are more easily com¬ 
pressed than sheared, whereas those with PR approach¬ 
ing 1/2 (large BjG) are rubber-like: they strongly resist 
compression in favor of shear. 

For most isotropic materials PR lies in the range 
0.2 < a < 0.5 Q. Materials with 0 < a < 0.2 are rare 
- Beryllium (a = 0.03), Diamond (a = 0.1) - and are 
very hard[l|. Typically, PR increases with pressure near 
linearly with the rate 6a/6P ^ 10“^ (GPa)“^ indicating 
a continuous loss of shear strength [JQ. 

An unusual pressure dependence of PR decreasing with 
rising pressure in solid hydrogen in the pressure range up 
to 24 GPa0 and solid helium up to 32 GPajsl was found 
by Zha et al. With the aim to investigate the distinc¬ 
tions in the response of quantum and classical solids to 
the applied pressure we calculated pressure dependencies 
of PR in the quantum (He, H 2 , Ne) and classical (Ar, Kr, 
Xe) cryocrystals under pressure. The calculations were 
performed using complementary semi-empirical (SE) and 
density functional theory (DET) with generalized gradi¬ 
ent approximation (GGA)approaches. The DET calcu¬ 
lations were performed using the EP-LMTO code RSPt, 
while the SE calculations were done using our own code. 
The calculation details have been published previously 
0- It is important to notice that the two approaches 
treat solid hydrogen in fundamentally different ways. SE 
approach deals with interaction between H 2 molecules, 
which are treated as nearly spherically symmetrical quan¬ 
tum rotators, while the DET can only treat fully ori¬ 
ented (classical) H 2 molecules, ignoring the zero-point 
rotations. The Pca2i oriented structure has been used 
for our calculations. 

One of the signatures of a quantum crystal is that 
it melts at temperature Tm much lower the Debye 
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temperature PH: 


QdITju > 1 . ( 2 ) 

Figure 1 shows pressure depeudeucies of 0^:) aud for 
solid hydrogeu, helium, ueou, aud argou. At zero pres¬ 
sure aud temperature the ratio Qd/T m is iufiuitely large 
for helium aud 8.5 for parahydrogeu. The ratio rather 
slowly decreases with risiug pressure. For example, at 
1 GPa it is still as high as 3.75 for helium aud 3.0 for 
parahydrogeu. For solid Ne at zero pressure ©d/T^ = 
2.7 thus makiug solid Ne a caudidate for the mauifes- 
tatiou of quautum effects. Other RGS, with ©^/Fm = 
1; 0.55; aud 0.35 for Ar, Kr aud Xe respectively, cau be 
regarded as esseutially classical solids. 

Rewritiug Eq. (1) iu terms of the ratio of the bulk (hy- 
drodyuami^i;^ to the trausverse (shear) souud velocity 
Vs we have [21 : 


{vb/vs)'^ - 2 
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The hydrodyuamic or bulk souud velocity vb cau be 
fouud from Equatiou of State: 


Vb = [dP/dp]^^‘^ = 
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(4) 


where P is pressure, p is molar mass, aud V is molar 
volume. Iu the calculatious of vb for H 2 we used our SE 
aud DET-GGA EOS from Ref. [l^, for He, Ar aud Xe 
from Refs. [li, [ill , aud for Kr from Ref. [Tol . We have 
also iucluded zero-poiut vibratious iu the Debye approx- 
imatiou iu our calculatious of p{V) aud vb{V) pH - 
Geuerally, to hud souud velocities vp aud vs oue has 
to hud a complete set of elastic moduli Cij. Iu the case 
of hep lattice there is a simplified scheme based ou lattice 
dyuamics[20|, HH, which makes it possible to circumveut 
the problem of calculatious of elastic moduli. Iu partic¬ 
ular, iu this approach it is possible to relate frequeucy u 
of the Ramau-active E 2 g phouou mode of hep lattice aud 
the shear elastic coustaut C 44 - 


C^i = ^--u\E2g), 

4V3aa 


(5) 


where a, c are the lattice parameters aud m is the molec¬ 
ular mass. The pressure depeudeucies of u{E 2 g) aud C 44 
were fouud for H 2 [13 aud hep RGS (hep He, Ar, Kr, aud 
Xe) [mil usiug both ab initio DET aud semi-empirical 
(SE) lattice dyuamics approaches. The shear velocity vs 
was obtaiued usiug the relatiou 


Vs = y C44/P, (6) 

where p is the deusity, disregardiug elastic the auisotropy 
of the crystal. A special case is solid Ne which preserves 
the fee structure up to at least 208 GPaji^ which makes 
the outliued procedure impossible. Eor this reasou for 
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FIG. 1. Debye temperature and melting temperature vs. 
pressure. Upper Figure: He and H 2 (insert); bottom Fig¬ 
ure: Ne and Ar (insert). Experimental melting curves: He 
and NeflH. Ar[l3]. H 2 melting curve corresponds to Kechin 
equation [Q. Debye temperatures were calculated using the 
many-body potentials: for He and H 2 pm , for Ne and Ar 

m. 


solid Ne we used results of lattice dyuamics calculatious 
by Gupta aud GoyaljUf. 

The souud velocities for H 2 aud He are giveu iu Ref. 
[H; the data for Ne by Gupta aud Goyal were published 
iu Ref. dil; the data for Ar, Kr, aud Xe will be published 
elsewhere. Pressure depeudeucies of Poissou’s ratios for 
helium, hydrogeu, ueou, argou, kryptou, aud xeuou cal¬ 
culated from souud velocities usiug Eq. (3) are showu iu 
Eigs. 2-5. 

Eigure 2 shows the pressure depeudeuce of Poissou’s 
ratio iu solid He obtaiued iu the framework of SE aud 
DET-GGA approaches iu comparisou with experimeutal 
results from Refs. Results which accouut for zero- 
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FIG. 2. (Color online). Poisson’s ratio of solid He as a func¬ 
tion of pressure. Theory: this work, Experiment: Zha et 
Nieto et al. m. 



FIG. 3. (Color online). Poisson’s ratio of solid H 2 as a func¬ 
tion of pressure. Theory: this work, Experiment: Zha et al.^ 


point vibrations (ZPV) and those obtained disregarding 
ZPV are presented. Both SE and DFT-GGA calcula¬ 
tions agree with the somewhat surprising experimental 
result of Poisson’s ratio decreasing with pressure. There 
is a reasonable fair agreement between the SE theoretical 
curve (comprising ZPV) and experimental data. Usually 
SE results are preferable at smaller pressures while at 
higher pressures the DFT approach works better. Com¬ 
paring the SE and DFT theoretical curves it is hard to say 
in which way the low-pressure SE results could continu¬ 
ously go over to the high-pressure DFT ones. It should 
be noted that the experimental points may show that 
around 30 GPa there is a minimum point at the pressure 
dependence of PR. 

As was said above, typically [El-Q Poisson’s ratio in¬ 
creases with pressure and tends to 1 /2 (the limit of zero 
compressibility) when pressure goes to infinity. It would 
appear reasonable to consider anomalous behavior of PR 
in such quantum solids as He and H 2 as a manifestation 
of quantum effects. Reasons for such understanding is 
the following. It is known that the He and H 2 lattices 
are swelled due to large zero-point vibrations (ZPV). If 
ZPV were not present, ’’classical” solid He and H 2 would 
have much smaller zero molar volumes (Uo^^(He) ^11.2 
cm^/mol; V T^(H 2 ) ~7.4 cm^/mol), i.e. the swelling ef¬ 
fect is huge [l8|. Until the volume reaches about the 
main effect of the external pressure is the suppression of 
the zero-point vibrations and not the compression of the 
electron shells. 

To check whether this explanation is correct we cal¬ 
culated PR of He disregarding ZPV, that is, for ’’classi¬ 
cal” He both in the SE and DFT approaches (dot-dash 
and dotted curves, respectively. Fig.2). As can be seen, 
the pressure dependence of PR with and without ZPV is 
qualitatively the same. Thus, the anomalous (descend¬ 
ing with rising pressure) behavior of PR is not a quantum 
effect. As can be seen from Fig. 2, the contribution of 
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FIG. 4. (Golor online) Poisson’s ratio of solid Ne as a function 
of pressure. Galculated using data on sound velocities by 
Gupta and Govalp^. 


ZPV into PR is positive. This fact is easily understood 
if we take into account that the introduction of ZPV is 
a step to liquation but PR of liquid is an upper bound 
for PR of any substance. Naturally, the relative value of 
this contribution increases with decreasing pressure and 
as pressure goes to zero it increases up to 15%. The ef¬ 
fect of ZPV is much higher in the case of ^He. Nieto et 
showed that the mixture ^He-^He has higher PR 
than pure ^He. For pure ^He they gave value of PR 0.473 
rather close to the liquid limit. 

The theoretical and experimental pressure dependen¬ 
cies of PR for solid H 2 are shown in Fig. 3. As can 
be seen, the SE and DFT-GGA approaches give the op¬ 
posite signs of the pressure effect on PR: PR decreases 
with rising pressure for SE and increases for DFT-GGA. 
Since in the experimentally studied pressure range (up 
to 24 GPa) the SE result agrees qualitatively with ex- 
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FIG. 5. (Color online) Poisson’s ratio of solid Ar, Kr, Xe as 
a function of pressure. Calculated using SE data on sound 
velocities (unpublished). 


periment [11, we conclude that at low pressures PR de¬ 
creases with rising pressure for solid H 2 . It is known that 
while the SE approach works well for molecular solids at 
low pressures, for higher pressures the DFT -GCA ap¬ 
proach is preferable. Thus the PR(p) curve can be sub¬ 
divided into three regions: At the low-pressure region SE 
is expected to work well, while at high pressures we can 
use the DFT-GGA approach. In the intermediate pres¬ 
sure range both approaches fail. The dot-dot dash curve 
shows schematically a possible continuous transition from 
the low-pressure asymptote to high-pressure one. Result¬ 
ing pressure dependence of PR for H 2 is non-monotonic 
displaying rather deep minimum. It should be noted that 
the transient region from the descending to the ascend¬ 
ing curves falls on phase II of the hydrogen phase dia¬ 
gram. As mentioned above, the SE and DFT approaches 
treat the orientational degrees of freedom in H 2 in com¬ 
pletely different ways: the former regards H 2 molecules 


as nearly spherically symmetric quantum rotators (as in 
phase I), while the latter considers classically oriented 
H 2 molecules (as in phase III), completely ignoring any 
quantum rotations or librations. It seems likely that this 
is the reason why SE and DFT give such drastically dif¬ 
ferent PR(p) curves for H 2 , while results for helium are 
qualitatively similar. It would mean that the PR mini¬ 
mum in hydrogen is related to the orientational transition 
at around 110 GPa, however more detailed study of this 
question is beyond the scope of the present work. 

A similar curve with a deep minimum was obtained for 
PR in solid neon (Fig. 4). The pressure dependence of 
PR was obtained from the SE theoretical results on sound 
velocities obtained by Gupta and Goyal[l9|. Unfortu¬ 
nately, experimental data on sound velocities in solid Ne 
exist for very narrow pressure range 5 - 7 GPaH^. In 
this region PR ^ 0.37. 

Figure 5 shows the pressure dependencies of PR ob¬ 
tained in the SE approach for Ar, Kr, and Xe. In con¬ 
trast with He, H 2 , and Ne, we obtained that PR for the 
heavy RGS increases with rising pressure. 

In conclusion, we present results of lattice dynamics 
calculations of Poisson’s ratio for solid hydrogen and rare 
gas solids (He, Ne, Ar, Kr and Xe) under pressure. Using 
two complementary approaches: lattice dynamics based 
on the semi-empirical many-body potentials and ab initio 
DFT-GGA we found three different types of the behav¬ 
ior of PR with pressure. While for solid He PR mono- 
tonically decreases with rising pressure, for Ar, Kr, and 
Xe it monotonically increases with pressure. For solid 
H 2 and Ne PR are non-monotonic with pressure display¬ 
ing rather deep minimums. To investigate the role of 
quantum effects we performed the calculations of PR dis¬ 
regarding zero-point vibrations and found qualitatively 
similar results, that is, we proved that the effects have 
a non-quantum origin, orblue We may rather say that 
the anomalies, discovered for H 2 , He and Ne, and quan¬ 
tum effects in these cryocrystals have common origins: 
weak intermolecular interactions and small masses of con¬ 
stituent atoms and molecules. 

We thank J. Peter Toennies for valuable discussions. 
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